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Abstract 

A relaxed notion of displacement convexity is defined and used to establish short 
time existence and uniqueness of Wasserstein gradient flows for higher order energy 
functionals. As an application, local and global well-posedness of different higher order 
degenerate non-linear evolution equations are derived. Examples include the thin-film 
equation and the quantum drift diffusion equation in one spatial variable. 
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1 Introduction 

In the last decade, the theory of gradient flows in the Wasserstein space has been a rapidly 
expanding area of research. With a wide range of applications to evolution equations and 
functional inequalities, this theory has received an extensive amount of interest. In this 
section we start by recalling some historical backgrounds of the theory and then we state 
the summary of our result. For a comprehensive discussion of all aspects of the theory we 
refer the reader to monographs [1] and |23j . 



1.1 Historical background 

The Wasserstein space V2{^'^) consists of the Borel probability measures on R"^ with finite 
second moment. The quadratic optimal transport distance, also known as the Wasserstein 
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distance 1^25 defines a distance function between any pair of measures /x, G V2{^"^) given 



where T{n, u) C P2(K'" x IR"*) is the space of probabihty measures with marginals /i and v. 
We will refer to such measures 7 as transport plans. 

It turns out that P2(IR™) has a rich geometric structure and a formal Riemannian cal- 
culus can be performed on this space. The first appearance of the Riemannian calculus 
on ^2(1^'") is due to Otto et al in [16] and [21]. It was shown in [21] that the solution 
of the porous medium equation dtu = Au^ can be reformulated as the gradient flow of 
the energy E{u) = J on the Wasserstein space. Since then, the interaction between 
the Riemannian space V2{^"^) as a geometric object and evolution equations as analytic ob- 
jects have attracted a lot of attention. This point of view is commonly called "Otto calculus". 

A notion which has been very important in the developement of this theory is the notion 
of displacement convexity. McCann in his thesis [19] introduced the notion of displace- 
ment convexity of an energy functional on the Wasserstein space. Under the displacement 
convexity assumption, he proved existence and uniqueness of minimizers of wide classes of 
energies, commonly referred to as potential, internal, and interactive energies. Displacement 
convexity had been defined before the development of the Wasserstein gradient flows, but 
after establishment of the Riemannian structure of the Wasserstein space, it turned out that 
displacement convexity can be interpreted as the standard convexity along the geodesies of 
the Wasserstein space. The displacement convexity condition, with its generalization to A- 
displacement convexity, has a central role in existence, uniqueness, and long-time behaviour 
of the gradient flow of an energy functional. 

Another important notion in the theory of Wasserstein gradient flows is the notion of min- 
imizing movement. Many of the rigorous proofs of the Wasserstein gradient flows are based 
on the method of minimizing movement. The minimizing movement scheme was suggested 
by De Georgi as a variational approximation of gradient flows in general metric spaces [12] . 
It was later used by Jordan, Kinderlehrer, Otto [16] and by Ambrosio, Savare, Gigli [1] to 
construct a systematic rigorous theory of Wasserstein gradient flows. This theory was soon 
used by many researchers to develop existence, uniqueness, stability, long time behaviour, 
and numerical approximation of evolution PDEs such as in [2], [5], [B], [S], [5], [TT], [13], and 



1.2 Summary of the results and outline of the paper 

In recent years, it has become apparent that Otto calculus also applies to higher-order 
evolution equations, at least on a formal level. The best-studied example is the thin-film 
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equation dtu = —V ■ (mVAm), which corresponds to the gradient flow of the Dirichlet energy 
E{u) = 1/ iVwpc/x. The hope is that gradient flow methods might help to resolve long- 
standing problems concerning well-posedness and long-time behaviour of this PDE. However, 
taking advantage of the gradient flow method has proved difficult. The main obstruction has 
been the lack of displacement convexity of the Dirichlet energy. The same problem arises for 
studying other energy functionals containing derivatives of the density. In [22] open prob- 
lem 5.17] Villani raised the question whether there is any useful example of a displacement 
convex functional that contains derivatives of the density. In [TU] , Carrillo and Slepcev an- 
swered this question by providing a class of displacement convex functionals. Therefore it 
was proved that there is no fundamental obstruction for existence of such energies. However 
because of the lack of displacement convexity, the Wasserstein gradient flow method has not 
been very successful in studying gradient flows of the Dirichlet energy and other interesting 
energies of higher order. 

Our result can be summarized as follows: 

• We introduce a relaxed notion of A-displacement convexity of an energy functional and 
in Theorem 12.41 we prove that, under this relaxed assumption, the general theory of 
well-posedness of Wasserstein gradient flows holds at least locally. 

• In Theorem 13. 6[ we prove that the Dirichlet energy, which is not A-displacement convex 
in the standard sense, satisfles the relaxed version of A-displacement convexity on 
positive measures. Hence the gradient flow of the Dirichlet energy is locally well-posed 
and the solution of the thin-fllm equation with positive initial data exists and is unique 
as long as positivity is preserved. 

• We show that the method developed to study thin-fllm equation applies to a range of 
PDEs of higher order and different forms. 

The paper is organized as follows. After recalling the backgrounds of the theory, in 
Section 2.2 we deflne the new version of A-displacement convexity which we call restricted 
X-convexity. Setting minor technicalities aside, the idea of restricted A-convexity can be 
summarized in two simple principles: Firstly, the modulus of convexity. A, can vary along 
the flow. Secondly, one can study A-convexity locally on sub-level sets of the energy. Note 
that the local analysis of gradient flows most likely fails without the help of energy dissipa- 
tion. For example, the Dirichlet energy is not even locally A-convex, because an arbitrarily 
small neighbourhood of a smooth positive measure contains measures with inflnite energy 
where A-convexity fails altogether. Instead, by taking advantage of the defining properties 
of the gradient flow, we study the flow on energy sub-level sets. The key observation is 
that typically flniteness of the energy implies some regularity on the measure which helps 
to elevate the formal calculations to rigorous proofs. For example, in 1-D, densities of flnite 
Dirichlet energy lie in H^. 

After deflning restricted A-convexity, we state our flrst result. Theorem 12. 4[ In this the- 
orem we prove that if an energy functional is restricted A-convex at a point /i, then the 
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corresponding gradient flow trajectory starting from exists and is unique at least for a 
short time. The proof is based on convergence of the minimizing movement scheme and the 
sub diff rent ial property that is carried over to the limiting curve. It is interesting that both 
of the constraints "locality" and "energy boundedness" are already encoded in the definition 
of the minimizing movement scheme (|T3|l . 

In Section 3 we apply the theory developed in the previous section to the Dirichlet en- 
ergy. We prove that the Dirichlet energy on M/Z, is restricted A-convex on the measures with 
positive density. This theorem re-derives the existing theory [5] of well-posedness of positive 
solutions of the thin-film equation by a direct geometric proof. To the best of our knowledge, 
this is the first well-posedness result for the thin-film equation based on Wasserstein gradient 
flows. Two key ideas are very useful in the proofs of this section: Firstly, the Wasserstein 
convergence and the uniform convergence are equivalent on energy sub-level sets. Secondly, 
finiteness of the energy can be used directly in the calculations of the second derivative of 
the energy along geodesies. 

In the final section, we show that the method developed in Sections 2 and 3 can be 
applied to a wide class of energies of different forms and of higher orders. Some important 
examples have been studied using this method such as equations of higher order of the form 
dtu = {—l)''dx{udl''~^^u), and equations of different forms, for instance the quantum drift 
diffusion equation dtU = —dxiudx^^^)- 

The Wasserstein gradient flow approach to PDEs has some interesting features. For 
example, it has a unified notion of solution which allows for very weak solutions and it is 
applicable to equations of higher order even with the lack of maximum principle. Also the 
minimizing movement scheme is a constructive method. Hence the proofs are constructive 
and one can derive numerical approximations based on the Wasserstein gradient flows similar 
to what has been done in [TT] and |13] . 



2 Well-posedness of the gradient flow 

In this section, we study the well-posedness problem of gradient flows on the Wasserstein 
space. Informally stated, a gradient flow evolves by the steepest descent of an energy func- 
tional. This idea can be formalized in several different ways, some of which carry over to 
general metric spaces. Here, we consider the Frechet sub different ial formulation of gradient 
flows. We will identify conditions on the energy functional that guarantee short-time exis- 
tence and uniqueness. The proof is based on a careful analysis of the minimizing movement 
scheme. 

Let us recall the notion of a gradient flow on a finite dimensional Riemannian manifold. 
The ingredients of a gradient flow consist of three parts: a smooth manifold M, a metric 
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and an energy E. Then the gradient flow of the energy E can be formulated as 

dfXt = Vt (velocity vector) 

Vt = —VE{xt) (steepest descent). 

Note that the role of the metric is to convert the co- vector dE = ^dx into the corresponding 
vector VE on the tangent space. 

In the case of Wasserstein gradient flows, the ingredients are given by: 7^2(1^™) as the 
manifold, the Wasserstein distance (and its infinitesimal version) as the metric, and an en- 
ergy functional as the energy. The formulation of a Wasserstein gradient flow is given in (fTOj) 
and as it can be seen, it is similar to its finite dimensional counterpart. 



2.1 Geometry of the Wasserstein space 

In this part we gather some basic elements of the Riemannian structure of the Wasserstein 
space P2(IR™)- Here, we work at a formal level and we refer the reader to [1] or [22] for 
rigorous proofs. It is worth mentioning that we consider the Euclidean space M™' as the 
underlying space of probability measures but one can replace M'" with any Hilbert space by 
slight modifications to the definitions as it is done in pQ. 

Consider the Wasserstein distance ([1]). The Brenier-McCann theorem |^ asserts that the 
minimum is always assumed and the minimal transport plan is concentrated on a graph of 
a map T;; : — > M™, provided that /i G (M™) where V^iW) is the set of absolutely 
continuous probability measures with respect to the Lebesgue measure. In this case, we have 
u = {Tl^)^fi, and one can rewrite the Wasserstein distance as 

W2ifi,u) = ([ |r;-/rffd/iy. (2) 

Assuming fi = udx and v = vdx are absolutely continuous measures, one can use the change 
of measure formula, given by the Monge— Ampere equation, to write an explicit relation 
between the densities u and v in terms of the optimal map: 

det{DTp{x) 

The optimal transport map also defines the geodesic fig between two measures /ig and fii 
given by the pushforward of the linear interpolation between the optimal map Tj^^ and the 
identity map: 

((l-s)/rf + .T^:;)^/io. (4) 

The appropriate class of curves inside ^2(1^™) which have a natural notion of tangent to 
them turns out to be the class of absolutely continuous curves ACfg^{[0, 00); 7^2(1^™))- 
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A curve fit belongs to ACf^^{[0, oo); P2(K™)) if there exist a locally L?'[dt) integrable function 
g such that 



W^2(/ia,/ife) ^ I 9{t)dt Va,6e[0,oo). 

J a 



The absolutely continuous curves are given by mass conservative deformations of the mea- 
sures i.e. they satisfy the continuity equation: 

9j/it + V.(/it^J =0 

for a velocity vector field Vj of deformations of fif This equation is assumed to hold in 
the distributional sense. In |1] Brenier and Benamou showed that V2{^"^) is a length space 
in the sense that the distance of two measures is given by the length of the shortest path 
between them: 



W2{fi,i^)=miU { \Vt\'dfit)-'dt s.t. dtfit + V.il^tVt) = 0; 1^0 = fi, fii = i^\ (5) 

where the infimum is taken over all curves in AC^([0, 1]; 7^2(1^™') )• For a given curve /it G 
AC^([0, 1]); "Pf (M™)) there might be many velocity vectors that satisfy the same continuity 
equation 9t/it + V.(/ifl^t) = 0. For instance vector fields of the form Ft + Vt where V .{ntFt) = 
all satisfy the same continuity equation. However there is a unique vector field that minimizes 
i.e. the one which defines the distance between fi and u. This optimal velocity vector 
field is defined to be the tangent vector field to the curve fit- The tangent vector field of 
fit can also be expressed in term of the optimal maps. If Vt is the tangent vector field of fit 
then 

T^*+^ - Id 

Vt = lim 



The converse is also true, i.e. for a given optimal map T^, the vector field T^ — Id is a tangent 
vector at fi for some curve that passes fi. The tangent vector fields are also useful in calcu- 
lating the derivative of the Wasserstein metric along curves. Let fit G /1C;^q^(M"'"; "Pf (M"^)). 
By [H Chapter 8] the derivative of the Wasserstein metric along the curve fit is given 
by 

^W^ifit. = 2 / {Vt, Id - rpdfit Wu G P2(M'") (6) 

at J^m 

where Vt is the tangent vector field to fit and (., .) is the standard inner product on M"^. 

The Wasserstein metric is closely related to a certain weak topology on 7^2(1^™); induced 
by narrow convergence: 

^„ =V fl ^ [ fdfir. ^ / fdfi V/ G , (7) 

JR'" JR™ 

where C^{W^) is the set of of continuous bounded real functions on M™. The topologies 
induced by the narrow convergence and the Wasserstein distance are equivalent for sequences 
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of measures with uniformly bounded second moments: 
lim 1^2 (/in, /i) = 



narrow 



jln > 

{fin} has uniformly bounded 2-moments. 



2.2 Wasserstein gradient flows 

Now we describe gradient flows on the Wasserstein space. Consider the energy functional 
E : V2{^"^) — [0, oo] and let its domain, D{E), be the set where E is flnite. Let /i lie in 
D{E) n V^iR""). A vector fleld ^ G L'^idfi) belongs to the subdifferential of E at if 

v&D{E) ^ ' 

We say that a curve [it G y4Cf^^(]R^, (^^'")) is a trajectory of the gradient flow for 
the energy i?, if there exists a velocity field Vt with |Vt|L2(^^^) e L;^^^(M~'~) such that 

dt\it + V ■ [[itVt) = (continuity equation), ^^^^ 
Vt G —dE{fit) (steepest descent) 

hold for almost every t > 0. We will refer to ffTOj) as the gradient flow equation. The con- 
tinuity equation, which is assume to hold in the distributional sense, links the curve with its 
velocity vector field and ensures that the mass is conserved. The steepest descent equation 
expresses that the gradient fiow evolves in the direction of maximal energy dissipation. 

Next we describe the link between Wasserstein gradient fiows and evolution PDEs. Let 
H = udx be in D{E) and let V G dE{fi) be a tangent vector field at fi. Consider a linear 
perturbation of fi given by the curve fi^ := {Id + eW)^fi for small values of e > where 
W G C^°°(M"*;M™). By the subdifferential inequality © we have 

E(u,)-E(u) f ,,,,,,, , ,. . ^E(u,)-E(u) 
lim sup —-^ — ^ / (V, W)udx ^ lim inf —. 

On the other hand, assuming regularity on and using standard first order variations 
we have 

E{u,)-E{u) f .5E{u)..du, 



lim = / ( \( ^udx 

where — r^-^ stands for standard first variations of E. Therefore 
bu 



[ {V,W)udx= [ ^^^d.ujx. 
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The continuity equation for the curve imphes that d^u^ = —V.^uW). Hence 

[ {V,W)udx=- [ {^^^V.{uW)}dx. 
Integrating by parts, we have 



Since W is arbitrary we have 

V(x) = V(^^^)(x) for/i-a.e. X. (11) 

OU 

Now assume that a curve /if = Utdx satisfies the gradient fiow equation (fTOj) . Steepest descent 
equation and (ITT]) imply 

V.i.) = -V(^)(.). 
By plugging Vt into continuity equation, we have 

a,„ = V, . (12) 

This is the corresponding PDE for the gradient fiow of the energy E. For example in the 
case of the Dirichlet energy E{u) = Jjg„ | Vwprfx, the first variation is given by ^ = -An. 
Therefore the corresponding PDE is the thin-film equation: 

dtu = -V.(mVAm). 

Our proofs are based on the minimizing movement scheme as a discrete-time ap- 
proximation of a gradient fiow which is described here. Let /iq € D{E), and fix the step size 
r > 0. Recursively define a sequence {M"}^^ by setting Mq = /iq, and for n ^ 1, 

M: = argmin (^(/i) + ^W^ (M^.^, /.) 1 . (13) 

The formal Euler-Lagrange equation for this minimization problem is given by 

U: G -dE {Ml) 

where t/^ = . Next we define a piecewise constant curve and a corresponding 

velocity field by 

/i[:=M;, Vr-=Ul, for (n-l)r<t<nr. 
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We have 



V; e -dE{fil) Vt > 0. 



(14) 



This equation suggests that /i[ is an approximation of the gradient flow trajectory of E 
starting from fio. 

There is a standard set of hypothesises that we assume throughout this section. We 
gather the hypothesises here: 

• E is nonnegative, and its sub-level sets are locally compact in the Wasserstein space. 

• E is lower semicontinuous under narrow convergence. 



The first two conditions guarantee existence and convergence of minimizing movement 
scheme ( lT3l) . The third condition ensures that measures of finite energy are absolutely con- 
tinuous, allowing us to use transport maps for studying the Wasserstein distance which 
simplifies the calculations, and allow us to view the sub differential as a tangent vector. Note 
that these conditions can be sharpened as in [1], but to make the presentation more appar- 
ent, we prefer to work in this more concrete setting. 

The main condition that guarantees existence and uniqueness of a Wasserstein gradient 
flow is given by the displacement convexity condition which asks for convexity of the energy 
along geodesies of the Wasserstein space. Let fig '■ [0, 1] — >■ V2(J^'") be the geodesic between 
/io, /ii € V2{X). An energy functional E is called displacement convex along fig if 



More generally the energy is called A-displacement convex or in short A-convex if the 
convexity is bounded from below by the constant A, i.e. 



Furthermore, assuming that E{fis) is smooth as a function of s, one can write a derivative 
version of A-convexity. In this case, E is A-convex along fig if 



E{fig) <: {1 - s)E{fio) + sE{fi^) se[0,l]. 



E{fig) <:{!- s)E{fio) + sE{^i^) - -s{l - s)Wi{^o. ^i) s e [0, 1]. 



(15) 




(16) 



An energy is called A-convex if it satisfies (fT5l) along all geodesies of V2{^"^)- 
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Figure 1: Restricted A-convexity 



2.3 Restricted A-convexity and local well-posedness 

Definition 2.1 (Restricted A-convexity.) We say that an energy E is restricted X- convex 
at ^ E D{E) with E{n) < c < +oo, if 36 > 0, such that E is X-convex along the geodesies 
connecting any pair of measures z/i, 1^2 G Bs{fi) fl Ec, where Bs{fi) = {v s.t. W^2(/^5 J^) < 6} 
and Ec = {v s.t. E^u) < c}. 

The following lemma has a key role in the arguments of Theorem 12.41 It shows how one 
can use restricted A-convexity assumption to study the subdifferential of an energy. 

Lemma 2.2 (Subdifferential and restricted A-convexity.) Assume that E is restricted 
X-convex at fi. Then a vector field ^ G L'^{dfi) belongs to the subdifferential of E at fi if and 
only if 

E{v) - E{^) ^ I (e, T; - Id) d/i + u) \/u e Bsif^) n Ec (17) 

where Bs{fi) fl Ec is the corresponding restricted X-convexity domain at fi. 

Proof. First we claim that for studying the subdifferential of the functional, it is enough to 
consider the restricted domain Bs{fi) fl Ec. Let ^ G L'^{dfi), we have to show that 



if and only if 



Ij, is trivial. 

For 1| assume that {z^n}]^ is a minimizing sequence for ( |T8|) . 
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In the case that hminf„^oo E{un) — E{fi) > we have 



W2{fi,l^n) \Jx ) W2ilJ.,l^n) 

E{U) - E{fi) 



+00. 



Therefore inequahty ( !T8|) is automatically true if liminf^^oo -E'(^n) — E{fi) > 0. Hence one 
only needs to consider sequences {i^n}'i' such that lim„_>.oo E{un) — E{fj.) ^ 0. Therefore, for 
large enough n we have E{vri) < c. On the other hand, z/„ /i. Hence f|T9|) and f lTSj) are 
equivalent. 

It is clear that ( ITTl) implies ( IT9l) . Conversely, let ^ G L'^{dfi) satisfy ( fT9l) . Let z/ G 
Bs{n) {lEc- Since £' is restricted A-convex at /i, we have A-convexity of E along the geodesic 
/is connecting /i to u. Therefore 

E{^s) ^ (1 - s)E{^) + sE{iy) - ^s{l - s)Wi{^l, v) Vs G [0, 1]. 
Dividing by s and reordering, we have 

^^^'^ ; ^^^^ ^ i^(^) - E{y.) -\{\- s)Wl{^^ V). (20) 
^ is in the sub differential of E at /i. Hence 



hminf ^^^^^ ^^^^ ^ lim 1 /" (e,rf -J^d/i 

5^0+ S s^0+ s Jx 

(e, t; - /d) d/i 



(21) 



where we used linearity of the interpolate map = Id + s{TI^ — Id). Therefore fl20|) and 
f l2T|) imply 

Eiu)-Eif^)^ J (e,T;-/d)d/i+ V|(/i,z/). 

□ 

The following lemma is used in Theorem 12.41 when we study weak convergence of tangent 
vector fields. 

Lemma 2.3 Let fx'l,i^it e AC^{[0,t\;V^{MJ^)) and let V^^ G L'^{dfx'l), Vt G L'^{dfXt)- Assume 
that 
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/ij fit uniformly on [0,t\. 

Vi'' weakly converges to Vt in the sense that \/U E C^{[0,t\ x M™), we have 



lim 

fc— >oo 





[V,',U{t,x))dfi'ldt= / / {Vt,U{t,x))dfitdt 



Then Wu G PslM"^) 

lim [ [ (V,\T\ - Id) dfi'ldt = [ [ {Vt,Tl-Id)dfjitdt. 

Proof. Since V^^ is weakly convergent by uniform boundedness principle we have 
supfc /o 4™ \ V^^\^dfi^^dt < +CX). Let 

M = sup /" ( /" \Vl'\^dfi'l+ [ \Vtfdfit)dt. 

k Jo JR™ jRrn 

Choose Tt G C°([0,f] x M™) such that 




\T'^^-Tt\^dfitdt<e\ 



We have 



Jo Jr™- \ * / Jo Jr"^ 

[ [ {Vi^Jd)di4dt- [ [ {Vt,Id)dfitdt 
Jo Jr'^ Jo Jr'^ 




+ 



r / {V,', T,> dfildt - f f {Vu T;,> dfjitdt 
Jo Jr"^ Jo Jr""- 



c 



We study each of the items separately. 

Since /i^ is uniformly converging to fit, the second moment of fi^ is uniformly bounded. 
In particular there is a compact set S* C [0,t| x R"* such that 



/if (it + sup / |x| fitdt) < e 

k Js'^ 



(22) 
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We have 



lim A = lim 

fc— >oo fc— >oo 



^ lim 



+ lim 

fc— >oo 



{Vt'',Id)dn^tdt- / / {VtJd)dntdt 

Jo JR™ 

Vt^Jd)d^4dt- I {Vt,Id)dfxtdt 
Js 

'yl',Id)d^Jldt- I {Vt,ld)dfxtdt 
Js'^ 



Because S is compact one can use weak convergence of V^^ on S. Hence the limit of the first 
term vanishes and we have 



lim A ^ lim 

k—>oo fc— >oo 



[Vt'^Jd) dfi'ldt 



(Vt, Id) djjLtdt 
Vtl^dfjitdt 

+ lim — ( I \x\'^dfj,'^dt + I \x\'^d^tdt 
k^oo 2e \Jsc Jsc 



^ - lim 

2 fc— >oo 



where we used Young's inequality with the constant e. By ( 1221) we have 

eM e 
hm A ^ h -. 

fc— >oo 2 2 

Since e is arbitrary we have limfc_j.oo A = 0. 



We now study B. Consider the measure 'y^ on 



given by 



It 



{idxT;,)#^4. 



Recall that the measure 7^ is the optimal plan with marginals /i^ and v. Since fi^ ^ fit, by 
the stability of optimal plans |23l Theorem 5.20], the set of optimal plans between fi^ and 
u is compact in the narrow topology and every limit point is an optimal plan between fit 
and u. On the other hand, because fit is an absolutely continuous measure, Brenier-McCann 
Theorem ensures that the optimal plan between fit and u is unique. This implies that the 
sequence 7^ converges narrowly for all t G [0,i\. Furthermore, the uniform convergence of 
fi^ implies that 'y^ have uniformly bounded second moment. We have 



{idxT;,))#fil 



narrow 



>^t = iIdxTn)^fit Vte[0,t]. 



(23) 



Taking the limit of B yields 



lim B = lim 

fc— >oo fc— >oo 



f [ lv,\T;, -Tt)dfi',dt 
Jo Jw^ ^ ' 



^ lim - 

fc— >oo 2 



eM 1 , 
^ 1 lim 

2 26 fe— >oo 



\Vl'\^dfi'idt + - 



\T;,-Tt\df4dt 



\T;,-Tt\^dfJidt. 
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by lifting to the optimal plans 7^^ = {Id x T^fe)^/i^, we have 



hm B ^ 1 hm / 

fc— ^00 

^ Jo ^R^xM"" 



y-Tt{x)fd^^dt. 



Since 7^ — 7* point-wise, 7^ has uniformly bounded second moment, and \y — Tt{x)\'^ is 
dominated by a constant times + + 1|, we can use dominated convergence theorem. 



Therefore 



lim 5 ^ — + — lim /" [ \y- Tt{x)\'^d-f^dt 

|y - Tt{x)\^d'-ftdt 
IT/" - Ttl^dfitdt 



eM 1 /•* 
eM 1 







2 2e 

eM e 

^ + -• 

2 2 

Finally, we study the last term C. We have 

{V,^Tt)dfi',dt 

{V,\T,)dfi1dt 



lim C = lim 



^ lim 

fc— >oo 



{Vt,T;:^)dfitdt 

{Vt,Tt) dfitdt 



+ [ [ {VuT;^-Tt)dindt . 
Jo Jr^ 

Since Tt G we can use weak convergence of V^^ for the first term. Hence 

Jo JR™ 



hm C ^ 

fe— >oo 



^ M 



T^^ - Tt\^d^tdt 



□ 



Theorem 2.4 (Existence and uniqueness of the flow) Let E : 7^2(1^'") — ^ [0, +00] 
6e a lower semi continuous energy functional with locally compact sub-level sets and let 
D{E) C V2{^"^)- Assume E{fi) < c < +00 and that E is restricted X-convex at fi. Then 
there exist t > and a curve fit G AC^ {[0,i\;V^(R"')) such that fit is the unique gradient 
flow of E starting from fi. 
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Proof. Let /i^ := be a piecewise constant solution to the minimizing movement scheme 
(IT^ with T = p The minimizing movement sequence is designed in a way that it converges 
to a hmiting curve in a very general setting. In [U Theorem 11.1.6] it has been proved that, 
under very weak assumptions which hold here, the minimizing movement scheme converges 
sub-sequentially to a limiting curve such that (after relabelling) Va > 

• /i*^ /it G AC^ ([0, a]; P2(K™)) uniformly in [0, a]. 

• The sequence {V^''} of the velocity tangent vectors to {/i^} converges weakly to Vt & 
L'^idfx) in Wx {0,T). 

• The continuity equation dtfit + V. (fitVt) = holds for the limiting curve. 

We need to prove that the limiting curve fit satisfies the steepest descent equation and that 
it is unique. Since /it is a continuous curve, we can find t such that /ij G Bs/4,{iJ,) for all 
t G [0,t| where 5 is the radius of restricted A-convexity at /i. We have 

E{fit) < limmiEijj'l) ^ E{fi) < c. (24) 

The first inequality follows from lower semi continuity of the energy and the second inequality 
follows from the structure of the minimizing movement scheme (fT3|l . Hence 

fit e Bs/iifi) n E, VtG[0,t]. (25) 

Since fit fit uniformly in [0,t|, we can find K G N such that W2{fit, fit) < ^/'^y ^ 
and Vt G [0,t|. Without loss of generality we assume that K = 1. Therefore ( 124|) and (!25|) 
imply 

The Euler-Lagrange equation of minimizing movement f|T^ implies that — G dE{fi^). 
Therefore using ( 126|) and the variational formulation of the subdifferential (Lemma 12. 2p we 
have 

E{u) - E{fi1) ^ [ {-Vt\ T;. - Id) dfil + ^WM. u) (27) 

for all V G Bs/2{fJi) H ^c- By construction, Bs/4{fit) H C Bs/2{fi'l) Ci E^. Therefore ([SZD 
holds for all u in Bs/4,{fit) H -Ec- By integrating fl271) over t and against a text function 
ij G C;?°((0,£); [0, oo)) we have 

f E{u)mdt- fE{fi^)mdt^ f I (-Vt\T\-Id)mdfi1dt + ^ / V|(/x^ i.)^(t)dt 
JO Jo Jo JM™ \ '^t / Z Jo 

(28) 

for all 1/ in Bs/i^fit) H -Ec- 

We take the limit of f l2S]) as A; — ?■ oo. By the lower semi-continuity of E 



fi'^ e Bs/2{fi) n E,. (26) 



E{v)i){t)dt 







E{fit)i){t)dt ^ I E{u)ip{t)dt- Mm ini E{fi^t)i){t)dt. 
Jo Jo 
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Lemma [2.31 implies that 

lim [ f (v^\T\-Id)tl,{t)d^Jldt= [ f {VuT;:^-Id)^{t)dfitdt. 

By the triangle inequality 

W.ifiuu) - W.ifitfit) ^ W^ifi^u) ^ W2{fit,i^) + W2{fi1,fit). (29) 

Therefore limfc_^oo W2{fif, v) = W2{fit, v)- Furthermore, since /ij uniformly, inequality 

implies that W2{fii, v) is uniformly bounded. Hence, by dominated convergence theorem 

lim ! Wi{iilv)mdt= [ W^{fit,u)mdt. 
Jo Jo 

In conclusion Wu G Bs/A{^it) H Ec and V-?/; G C^((0, i); [0, oo)) we have 

f E{u)mdt - f E{fit)mdt ^ f I {-Vu - Id) mdiitdt +^ f W^i^it, u)i;{t)dt. 

Jo Jo Jo JM™ ^ Jo 

(30) 

Let to be a Lebesgue point of the map t h-)- E{^t)ip{t)dt + (Vt^T^^ — Id)dfXt + 

f Io^2il^t,i^)ipi't)dt. By considering a sequence of smooth moUifiers ipn converging to the 
delta function at to; the inequality ( 130|) is reduced to 



if^to,^)- (31) 



Therefore ( 13T]) holds for almost all t. By Lemma [2.21 we have 

Vt G -dE{j2t) for almost all t G [0,t]. 

We now study uniqueness of the solution. The available uniqueness proofs in the case of 
A-convexity can be repeated in the domain of restricted A-convexity because as soon as the 
flow exists clearly it dissipates the energy and the trajectories of the flow starting from /i re- 
main in the domain of restricted A-convexity for a short time where one can use A-convexity. 
Hence uniqueness arguments can be repeated similar to the available proofs such as in [H 
Theorem 11.1.4]. Therefore we provide only the key ideas here. 

Assume that we have two gradient flows fi} and fif both starting from fi. One can show 
that W2{fJ^l,fit) is absolutely continuous in time and 

dt h ™ h 

By differentiability of Wasserstein metric ([6]), for almost all t G [0,t| we have 

2 h^o h Lm \ * 
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Therefore 

/i?) ^ £^ (V,\ Id - T;|) rf/i? + (V/, Jrf - T^l) (32) 
Consider (13T|) along /xj. We have 

E(/i?) - E{^^\) ^ j^^ {-Vt,Tll - Id) df^l + V2(/x,\/x?). (33) 
Rewriting (133|) again along /i^ and using (!32l) result in 

Hence 



3 Wasserstein gradient flow of the Dirichlet energy 

In this section we prove a local well-posedness result for the gradient flow of the Dirichlet 
energy on S^. Energy or in short E in this section always refers to the Dirichlet energy 

i -(d^uf if fi = udx,ueH\S^), , , 

E(/i) = <^ 2^ ^ ^ ^ ' ^ (34) 

[^+00 else. 

When convenient, we refer to an absolutely continuous measure /i = udx by its density u. In 
particular, by a smooth or positive measure we mean a measure with a smooth or positive 
density. 

3.1 Dirichlet energy on 

The underlying space of the measures that we study in this section is . We identify 
with M/Z. Because is a manifold, the theory developed in the previous section should be 
slightly modified. On a Riemannian manifold, there is the issue of existence and regularity 
of the optimal maps. This question has been an active area of research. Ma-Trudinger-Wang 
condition in [18] is a famous example that studies this issue. In the case of S**", the problem 
has been addressed and positive results are available such as |T5] and [17] . The results guar- 
antee that between any pair of smooth positive measures /i, v on S*"* there exists a unique 
smooth optimal map T^. To apply Theorem 12.41 to 5*^, one has to replace the inner product 
in the subdifferential definition ([9]) with the inner product on the tangent space of . For a 
given optimal map T, the distance |T(x) — x| might not coincide with the geodesic distance of 
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S^. As was suggested in [10], this problem can be solved by representing T : [0, 1] — t- [— |, |] 
where T(x) is the smallest element of T(x) G M/Z such that |T(x) — a;| ^ | and by rela- 
belling by [T(0),T(0) + 1]. It is easy to check that T is an optimal map from [0, 1] to 
[T(0),T(0) + 1], furthermore T is monotone and the geodesic distance of coincides with 
|T(x) — x\. 

Let fiQ = u^dx and yUi = Uidx be two smooth measures on and let T be the optimal 
map between them. Then by Monge- Ampere equation ([3]) we have 

Ui{T[X)) 



T'(x) 



Since the geodesies are given by the push forward of linear interpolation of the optimal map 
and the identity map, the explicit form of the geodesic Ug between uq and Ui is given by 

In the notation of the previous section, if we think of / = T — Id as the tangent vector field 
that connects uq to Ui, the geodesic equation can be written as 

usix + sfix)) = -^^^ (35) 
1 + sf'{x) 

We will see in Lemma 13.31 that for studying restricted A-convexity of the energy, it is 
enough to consider measures with smooth and positive densities. By the derivative formu- 
lation of A-convexity f|T6l) the energy is A-convex along the geodesic fis connecting fiQ to fii 
if 

^E{u,)^XW^iuo,u,). (36) 

We start by proving that the Dirichlet energy is not A-convex for any A. This is known 
to the community, and in |T0] Carrillo and Slepcev proved that the Dirichlet energy is not 
convex on . We will study a scalable family of functions to prove the lack of convexity for 
any A G M. Let Us be the geodesic connecting uq = u to ui, and let / be the corresponding 
tangent vector field. We compute the second derivative of the energy along the geodesic 



_ d^ 






/ 

s=0 JS 



{dyUs{y)fdy. 



di' 

By the Monge- Ampere equation fl35l) and the change of variables y = x + sf{x) we have 
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s=0 



ds^ 
ds^ 
ds^ 



dx d u{x) 
oJs^ \dydxl + sf'{x] 



dy_ 

dx 



dx 



uix] 



Js^ \l + sf'{x) n + 

r {{l + sf '{x))u'{x) - su{x)r{x)f 
oJs^ 



(1 + sf'{x)y 
{{f"uf + S{f"u){f'u') + 6{f'ur}dx 



1 + sf'{x)) dx 
dx 



(37) 



If the energy E is A-convex, then 



{{f"uf + 8iru)ifu') + 6{fu')'}dx^X [ fudx. 
1 Js^ 



(38) 



B 



We use fl55]) as a guide to find a counter-example. In the example, u and / are not smooth, 
and we cannot directly apply this computation but. Lemma [3.31 validates the calculations. 

We view as the interval [—1/2, 1/2] with the endpoints identified. The construction 
of the example is simple: let m = 1 — 4|x| and /' = u~^, this forces the integrand of A to be 
negative, and the rest follows from a scaling argument. We have to make some modifications 
to the functions so that the integral converges and the mass is normalized to 1. We define u 
and /' as follows 



u{x) 




^ Ixl ^ 



^ \x\ ^ 
^ \x\ ^ 




^ kl ^ 



^ \x\ ^ 
^ 1x1 ^ 




Figure 2: Graph of u and /' 
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By scaling Uh{x) := hu{hx) and /^(x) := j^f'{hx) we have 

A = -Cih^ , B = C2 

for some positive constants Ci,C2- If the energy is A-convex then we must have A ^ XB 
and it should hold uniformly for any such u and /. But for a fixed A, we can choose h large 
enough so that the opposite inequality holds. This means that the Dirichlet energy is not 
A-convex on V2{S^)- 

In the example above, by pushing h to larger numbers the lack of convexity becomes 
worse. By looking at the equation of u{x), it is clear that the Dirichlet energy of u{x) gets 
bigger for larger values of h. This example hints that one of the obstructions against the 
A-convexity of the Dirichlet energy is the magnitude of the energy which can be controlled 
on energy sub-level sets. 



3.2 Restricted A-convexity of Dirichlet energy 

Lemma 3.1 (Uniform convergence on energy sub-level sets.) Uniform convergence and 
Wasserstein convergence are equivalent on energy sub-level sets of Dirichlet energy on . 
In particular for two measures fii = Uidx and fi2 = U2dx with E{fii), E{fi2) < c < +00 we 
have 

W^{lJul^2)>a\ui-U2\i (39) 
where a = a{c) and /3 are constants. 

Proof. One side of the equivalence is easy. Assuming Un "™'^°'''") mq we have 

ipUndx — > / il'udx G C'^{S^) 

which implies Wasserstein convergence of /i^ = Undx to fi = udx by ([H]) and finiteness of the 
second moments on S"^. 

For the converse inequality, we first study the regularity of a measure with finite energy. 
Let V = vdx G Ef.. By Poincare's inequality and J^^ vdx = 1 we have 

\v\'^dx ^ / \v'\'^dx + 2 / \v\dx+ / dx^c + 3. 

Therefore if^-norm of v is bounded by its energy. The Sobolev embedding theorem implies 
that V is (7°'^/^ continuous and we have 

1/2 / N 1/2 



v{x)-v{y)\^i \v'\'^dx\ i dxj ^ ^/c\x-y\ \/x,y e S\ (40) 
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Figure 3: Wasserstein 



Uniform 



Therefore the modulus of continuity is a/c. Let fii and fi2 be as in the assumption. Therefore 
Ui and U2 are C"'^/^ continuous with constant y/c. Assume that \ui — U2I00 ^ h > 0. In 
particular without loss of generality assume that for some point Xq G we have Ui{xo) — 
^2(2^0) ^ h. For every x G S*^, we have 



Ui{x) ^ ui{xo) — \/c\x\ 

U2{x) ^ U2{xq) + a/c|x| 



(41) 



Therefore Ui lies above and U2 is below the star-like shape in Figure (|3]). Call the star-like 
shape by S. Consider a rectangle R in the center of S with height | and width | where k is 
the width of S at the height U2(xo) + |. We have k = ^ and the area of R is given by 
In order to transport the measure fii to fi2, some mass at least equal to the area of R should 
be transported outside of S. Therefore 



H^2^(/xi, /i2) ^ {area of R}{distance required to move R outside of S}' 



48c'V 
1 

384? 



\Ui - U2 



7 

00 ■ 



□ 



Lemma 3.2 (Lower Semi Continuity.) The Dirichlet energy is lower semi continuous 
with respect to the Wasserstein metric on S^. 
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W2 



Proof. Let u„ u. Since we have to prove E{u) ^ liminf„_j.oo -E'(un), we can assume 
that E{un) is bounded. This imphes that if^-norm of the sequence is bounded. By Ba- 
nach-Alaoglu theorem, Un has a weak hmit point v G H^. Therefore u„ converges to v 
strongly in L^. We claim that v = u. Let if) G C^{S^) we have 



il){x){u — v){x)dx 



1p{x){u — Un){x)dx 
1p{x){u — Un){x)dx 



1p{x){Un — v){x)dx 
+ (/ \i,\Hxf'\! \Un-v\^dxf'^ 



By Lemma 13.11 the Wasserstein and uniform convergences are equivalent on energy sub- level 
sets. Therefore the first term in the last inequality goes to zero. The second term also goes 
to zero because m„ converges to v strongly in L^. Hence u = v almost everywhere. Because 
u and V are continuous, we have u = v. The Dirichlet energy is known to be lower semi 
continuous under weak convergence (for example see p3l Theorem 8.2.1]). Hence we 
have Eiu) = E(v) ^ liminf„ Eiur,). □ 



The following lemma validates smooth calculation in the sense that for studying restricted 
A-convexity of the energy, one can study restricted A-convexity of the energy only on smooth 
measures. 



Lemma 3.3 (Approximation by smooth measures.) Let fi G Ec- Assume that the en- 
ergy is restricted X-convex on smooth measures in Bs HEc- Then E is restricted X-convex at 
/i. 

Proof. Let /xq = uodx, /ii = uidx G Bs H Ec and let rjk be a standard smooth moUifier 
converging to the Dirac delta function. Define Uk,i{x) := rjk* Ui{x) for i = 0,1 where * is the 

convolution on S^. Since Uf^^i ""''^°'''"'^> Ui, by Lemma [3. II we have Uk,i — ^ Ui. Therefore for 
large enough k we have Uk^t G Bs{u). The energy also converges, because 

E{uk,i) = / {d^{ui * r]k)f dx = / {{d^^Ui) * rjkf dx > / (d^Uif dx = E{ui). (42) 

Hence Uk.i G Bs{u) (lEc for large enough k. By smoothness of Uk.i and the assumption of the 
lemma, we have A-convexity of the energy along the geodesies Uk^s connecting Ukfl to Uk,i 

E{uk,s) ^ (1 - s)E{uk,o) + sE{uk,i) - ^s(l - s)W^{ur,,o, Uk,i). (43) 

Let 7fc and 7 be in order the optimal plan connecting /x^ = Ukfldx to ^k,i = Uk^idx and the 
optimal plan connecting /ig to /ii. By stability of the optimal plans [231 Theorem 5.20] 7^ 
converges in narrow topology to 7 along a subsequence which after relabelling we assume to 
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be the whole sequence. Equivalence of narrow and Wasserstein convergence (l39l) on x 
implies 

/ifc,, = ((1 - s)ni + sn^)^ 7, ((1 - s)ni + sn^)^ 7 = /i, 

where 11* is the projection to the i*'* coordinate and fig is the geodesic connecting fiQ to fii. 
The lower semi-continuity of Dirichlet energy 1X2) yields E{fis) ^ liminffc_^+oo -E'(yLtfc,s)- Hence 
by taking the limit of ( H3l) we have 

E{ix,) ^ (1 - s)E{fXo) + - ^.3(1 - s)W^|(/io,/ii). 

□ 

In the following lemma we prove that the energy is finite along a geodesic, provided that 
the energies of the end points are finite. 



Lemma 3.4 (Energy of the interpolant.) Let fiQ = u^dx and fii = udx he two smooth 
measures with E{fio), E{fii) < c < +oo and Uq,Ui > m > 0. Then there are constants 
c < +00 and m > depending only on c and m such that E{fis) < c and Ug > m along the 
geodesic fig connecting fiQ to fii. 

Proof. By C*°'^/^ continuity of the densities, there exists M = M(c) such that uo{x),ui{x) < 
M for all X G S*^. Let T : — > be the the optimal transport map between uq and Ui. 
By Monge- Ampere equation ([3]) we have 



\T'(x)\ ^ 



M 



m 



By taking the derivative of Monge-Ampere equation ([3]) we have 



\T"(x) 



u'q{x)ui{T{x)) - uo{x)T'{x)u[{T{x)) 



ui{xy 



M 



(44) 



Now let fig be the geodesic connecting /io to fii. We have 

uo{x) 



Us{{l - s)x + sT{x)) = 
Plugging in bounds on T' yields 

Us{x) ^ 



;i - s) + sT'{x)' 



(45) 



m 



Hence Us > rh where m = rh{c,m). Taking derivative of the equation (l45ll we have 

u'q{x)[{1 - s) + sT'{x)] - suo{x)T"{x) 



Ug{{l — s)x + sT{x)) 



((1 



sT'{x)Y 
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Using the bounds on T' and (jH]) we have: 



;i - s) + sT'(x) |2 I (1 - s) + sT'{x) |3 

^ (-)Vo(^)i + {-)'\uo{x)\ + {-rK{x)\. 

m m m 

Taking integral from both sides yields E{us) < c where c = c(c, m). □ 

The idea of the following lemma was suggested by my supervisor Almut Burchard. This 
lemma will be used in calculations of the second derivative of the energy in Theorem 13.11 

Lemma 3.5 (Interpolation inequality.) For every a G M there exists a constant A < 
such that 

\f"\L^ - c^irioo -Mf\L^>o V/ e C^iS'). 

Proof. Consider the Fourier expansion f{x) = Xlfcez '^fc^*^^'^^- have 

l/'loo = sup I ^i2A;7rafce^2'^^^'| ^27r^|A;afe| = 27i^\k^al\^\al\^\h^ . 



^^•^^ kez kez 



Holder's inequality with exponents jq, and ^ yields 



vfcez / \kez / Vfcez 



The term 27i{J2kez 1^1 ^•'^^^ = d is a constant independent of Cfc. Therefore 



I/loo ^ d\fX!'\f\y' 



lL2 \-l lL2 • 

By the arithmetic-geometric inequality for a constant /3 we have 

\f'\oo^d\f"\%'\f\'l!' 



4 dci(5~^ 
Putting 13 = {^ad)~ and A = yields 



5 

"irioo^iru^-Ai/u^. 



□ 



We are now ready to prove the main theorem of this section which shows that the Dirichlet 
energy is restricted A-convex at positive measures. 
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Theorem 3.6 (restricted A-convexity of the Dirichlet energy.) Let fi = udx he a mea- 
sure with E{u) < c < +00 and u > m > 0. Then 3A = Xc^m such that E is restricted X-convex 
at fj,. 

Proof. We first claim that the second derivative of the energy at a positive measure v is 
uniformly bounded from below along any smooth vector field. Let v = vdx be a measure 
with E{v) < c and v > m. Let ui be another smooth measure and let / be the vector field 
defining the geodesic Ug = {Id + sf)^v that connects v to vi. By ( 1371) we have 



d'^Eiu, 



ds^ 



= 2 / {f"vy + 8{vf"){v'f') + Q{v'frdx. 

s=0 JS^ 



Recall that W|(z/, ui) = f^^ f{xYv{x)dx. By fl36|) the energy is A-convex at v, if for all such 
vector fields 

2/ {vf"f + ^{vf"){v'f) + Q{v'ffdx-xj vfdx^O. 



Jsi 

By completing the squares we have 

2/ {ivf"f + 8{vf"){v'f') + 6{v'f'f}dx-x[ vfdx^ f {/"V - 52/'V2}dx - A / vfdx. 
The lower bound on the density v > m yields 

/ {f"\'' - 52f%'^}dx - A / vfdx ^m^ [ f'^dx - 52 / fVdx - mX [ fdx. 
Holder's inequality and energy bound E{u) < c imply 



I f'^dx- 52 [ fu'^dx-mX [ fdx^m^ [ f'^dx - 52c\f'^\oo - mX [ fdx. 



By reordering and absorbing the constants in A, the energy is A-convex along Us at u if 
V/ G C^{S^) we have 

irii^-«l/'lL-A|/|i.^O (46) 
where a = By Lemma [3.51 the claim has been proved. 

Now consider the energy sub-level set E^. By Theorem 13.11 Wasser stein convergence im- 
plies uniform convergence on Ec. Therefore there exists a. 6 = 6c such that we have v > m 
for all z/ = vdx E EcH Bs{fi). Assume that I'o^'^i & E^H Bi^v). Let be the geodesic 
connecting to Vi. By Lemma [3.41 there exist rh and c depending only on c and m such 
that E{vs) < c < +00 and Vg > m > 0. By the argument at the beginning of the proof there 
exists a A = Am,c such that E is A-convex along the geodesic Ug. The constant A is uniform 
for all pairs of smooth measures inside Ec fl Bs{fi). Therefore, by Lemma [3.31 is restricted 
A-convex at fi. □ 
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Corollary 3.7 The gradient flow trajectory of the Dirichlet energy on with a positive 
initial data exists and is unique at least for a short period of time. 

Corollary 3.8 The positive periodic solutions of the thin-film equation dfU = —dxiud^u) 
are locally well-posed. 

4 Other classes of equations 

In this section, we show that the theory developed in the last two sections can be applied 
to a wide class of energy functionals and evolution equations of higher order and different 
forms. Note that the result of Theorem 12.41 is general and it can be applied to any energy 
functional, provided that it is restricted A-convex. The corresponding lemmas from Section 
3 for the energies studied here can be derived in a similar fashion with minor modifications. 
Hence, we discuss the proofs only briefly 



4.1 Higher order equations 

The family that we study here is of the form E{u) = ^ J^i \u^''^\'^dx for k E N. The flow of 
this family of energies corresponds to the solution of the higher order non-linear equations 
of the form dtu = {—!)'' dx{udl''~^^u). 

Consider u G D{E). Finiteness of \u\ijk in particular implies that the H^-noim. of u 
is bounded. Since we only used the if^-norm bounds in Lemmas 13.11 13.21 and 13.31 they 
automatically follow for this class of energies. Therefore, Wasserstein and uniform conver- 
gence are equivalent on energy sub-level sets, E is lower semi continuous, and one can use 
approximation by smooth functions to study convexity. 

In Lemma 133] we derived bounds on T" by taking derivatives of the explicit formula of T' 
given by the Monge-Ampere equation. In the same fashion, one can find bounds on higher 
derivatives of the optimal map by taking more derivatives of the Monge-Ampere equation. 
For generalization of Lemma 13. 5[ we have to show that Va 3A such that 

|/(-+i)|2^_«|/M|^_A|/|2, ^0 

for every smooth vector field /. By induction assume that for any am > there exists 
Am ^ such that 

I/^'^^IL^— l/^'^+'^li^-A^I/li.. (47) 
Let am+i > be given. By applying Lemma [3.51 to there exists a A ^ such that 

|^C^+i)|2^ ^ ^If^'-^'^ - A|/(-)|i.. (48) 
2am+i 
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Put am = — 2A, by (H7|) there exists ^ such that 

2a 

|/M|2^ ^ on M/Z. Therefore 

Plugging into ( HHl) yields 



2am+i 2 

Therefore 

l/^'-^^^lL-^l/^'-^^^li^ ^ ^l/^"^^^^li^ + AA„|/ii.. 

2 20Lm+\ 

By ^ and by setting A^+i = ^AA^, we have 

Va^+i > 3A„+i ^ s.t. |/(-+i)|^^^|/(™+2)|2^_^^^^|y|2^_ (49) 

In conclusion, all the Lemmas in the previous section can be applied to higher order 
energies. We now study convexity of the energies along smooth vector fields on a measure 
\x = udx with positive density u > m and finite energy E{u) < c < oo. 



ds'^ ^ ^jgi^'^dydx' l + sf'{x)^ dx 

Since we study all the different orders at the same time, we consider the general form given 
by a polynomial P which is determined by the order of the energy. We have 

^,\s=oE{u,)= I |n/(^)p + P(n,nW,...,n('^);/,/«,...,/(^-i))cix 

where P is of order at most 2 with respect to each of its entries, and the order of the 
derivative of each term in P is at most k. At a measure with positive density and finite 
energy, we have u > m and |m^*^|oo < M for all i < k where M depends only on E{u). Also 
^ |/^'^~^''|oo for alH < A; — 1. Therefore similar to the calculation of Theorem 12. 4[ for 
positive constants /3i,/32 we have 



2 



d_ 

ds' 



oo ■ 
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Therefore E is convex at u if we can find A such that 



L2 

This imphes that the energy is A-convex at u because by poj) for a = am,c there exists A 
suicti ijliciij 

|/W|2^ - a\f^'^-'^\l - A|/|i. ^ V/ G C^{S') 
Hence we have proved the following theorem. 

Theorem 4.1 The energies of the form 

/gi \dl;u{x)\'^dx /i = udx, u G H'^{S^), 



E(u) 



-oo else. 



are restricted \-convex on the positive measures with finite energy. In particular, periodic 
gradient flow solutions of 

with positive initial data exist and are unique for a short time. 
4.2 Different forms of equations 

Consider the energies of the form E{u) = J^j g{u, d^u)dx. We start by calculating the second 
derivative of the energy along a geodesic induced by a vector field / G C°°{S^). 



ds 



2 ^s=oE{Us 



d' 






s=0 


d^ 






s=0 


d^ 






s=0 



u{x) dx d u{x) 



1 \1 + sf'{x) ' dy dx 1 + sf'{x) J dx 
u(x) 1 _ , u(x 



dy , 
—dx 



Is." \l + sf'{xyi + sf'{x) ^'l + sf'{x) 
where we used the change of variable y = x + sf{x). Therefore we have 



;i + sf\x))dx 



^.UM^Jjf f"]A 



f 
f" 



dx 



(50) 



where the matrix A is given by 
A = 



^2^(0,2) 



Note that if A is positive definite, then the energy is convex. We study the class of 
the form g{u,dxu) = \dx{u°')\'^ with a > 0. Finiteness of the energy implies that u"" is 
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C^'^ continuous with modulus of continuity smaller than the energy. Because a > 0, m is 
continuous and since f^^udx = 1, there exists a point Xq with u{xo) = 1. Without loss of 
generality we assume Xq = 0. We have 

\u{xy - 1| ^ \/c|a?[ =^ u{x) ^ (1 + \/c\x\)^. 

Therefore there exists a uniform M < oo such that u < M for all u G Ec- We now briefly 
discuss the corresponding lemmas from Section 3. 

Equivalence of Wasserstein and uniform convergence on energy sub-level sets. 

Let U2{xq) — Mi(xo) > h. Then we have U2{x) ^ ^2(^0) ~ a/c|x — xo| and ^ 
u\{xq) + \/ c\x — xo|. Therefore the star-like shape in Lemma 13.11 should be replaced by 
a modified version, given by {u2{xq) — ^/d\x^^x^\)'^ and {u1{xq) + \/ c\x — xo|) and the 
rest of the proof goes similarly. Hence, we have equivalence of the Wasserstein and uniform 
convergence on the energy sub-level sets. 

Lower semi-continuity and smooth approximation. Having Lemma 13.11 for this 
class of energies, the proof of Lemmas 13.21 and 13.31 can be repeated by replacing u with u". 
Hence the energy E{u) = J^^ \dxU°'\dx is lower semi continuous and one can use approxima- 
tion by smooth functions. 



Energy of the interpolant. Let u be bounded away from zero m > m > 0. When 
a ^ 1 

m'^""'^ / \d,u\''dx ^ E{u)= I m2(«-i)|5^^|2^^ ^ m2(«-i) / \d^u\''dx (51) 



'SI 

and when < a < 1 



M2(»-1) / \d^u\^dx ^ E{U)= [ M2(«-l)|5^.^|2^a; ^ ^2(a-l) f \Q^^\2^^^ (52) 

By equivalence of Wasserstein and uniform convergence, there exists 6 such that v > m for 
all V G Bs{u) n Ec- Also we have proved that v < M for all v G E^. Therefore we can 
refer to the calculation for the Dirichlet energy and just compare the energy of the geodesic 
with the corresponding Dirichlet energy using (15T|) and ( l52l) to find a bound on the energy 
of interpolate points along a geodesic. 

In conclusion, all of the required lemmas are true. By (ISOll . along a geodesic induced by 
a smooth vector field / we have 

^l^,E{us) = 2aV('^-i) {{uf"Y + 4(1 + a){un{u'n + (1 + a)(l + 2a){u'f'f) dx 
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for some constants ai,a2. Similar to ( H6|) . we have 

^l^,E{u,)-m\ f ufdx ^ \f"\l,-a\r\l^-\\f\l. ^ 

where the last inequality follows from Lemma 13.51 We have proved the following theorem. 
Theorem 4.2 For every a > 



E{u) 



-oo else. 



is restricted X-convex on positive measures with finite energy. In particular, periodic gradient 
flow solutions of 

dtu = ~2ad^{ud^{u''-^dlu'')) 
with positive initial data exist and are unique for a short time. 

An interesting example is the Fisher Information 

Eiu) = ]- f \d,,u{x)^^dx 
which corresponds to the quantum drift diffusion Equation 

^2 

dtu = -d^{ud^ ^ 



u 



Therefore we have local well-posedness of periodic solutions of the quantum drift diffusion 
equation with positive initial data. 

Another interesting case is the limiting case a = 0. The corresponding energy can be 
written as E{u) = | J^i \dxlogu\'^dx. Finiteness of the energy result in C°'2 continuity of 
logM. All of the lemmas can be repeated in a similar fashion for this energy. Furthermore, 
finiteness of the energy implies a lower bound for the measure because 

|logu(x) -log(l)| ^ =^ e"^^n(x)^e^. (53) 

Therefore positivity is preserved along the flow. By ( 150|) we have 

^,\s=oE{fis) = [ 2{fr + 4{f"){-f')-2{-frdx. 
ds"^ Jgi u u 

By 13.51 there exists A such that 

^^s=oE{^^s) > irii^-eceV^I/'lL-AI/U^ ^ 0. 



ds 

Hence E is restricted A convex at n G E^.. Furthermore, since there is a uniform lower bound 
g-v^ all V E Ec, the constant A is uniformly bounded along the flow. Therefore, the 
gradient flow is globally well-posed and we have the following theorem. 
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Theorem 4.3 Wasserstein gradient flow of the energy 

^ \^Jgi\dxlogu\'^dx fi = udx, u E H\S^) 
1 +00 else. 

is globally well-posed. Hence the equation 

dtu = -d^{udl-^) 
with periodic boundary condition is well-posed. 

Remarks. There are some simple and some more challenging directions to extend the 
developed method to other classes of equations. As a simple application, one can construct 
other classes of restricted A-convex functionals by combining the ones already studied. For 
example, the solution of the energy E{u) = Jgi{\dxu\'^ -\-e^}dx, which is the Dirichlet energy 
with a perturbation, is globally well-posed. The reason is that the second term forces the 
energy to remain positive. One interesting problem is the analysis of equations in higher 
dimensions. Our method is utilizing Sobolev embedding theorem on energy sub-level sets 
which is getting weaker on higher dimensions. An interesting question is whether it is pos- 
sible to solve this problem with studying higher order energies. 
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